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ABSTRACT 

We present the results from a survey of z-dropout objects selected from ~ 1550 
deg^ of multicolor imaging data from the Sloan Digital Sky Survey (SDSS), to search 
for luminous quasars at z 5.8. Objects with i* ~ z* > 2.2 and z* < 20.2 are 
selected, and follow-up J band photometry is used to separate L and T type cool 
dwarfs from high-redshift quasars. We describe the discovery of three new quasars, 
SDSSp J083643.85+005453.3 {z = 5.82), SDSSp J130608.26+035626.3 {z = 5.99) and 
SDSSp J103027.10+052455.0 (z = 6.28). The quasar SDSSp J083643.85+005453.3 is 
a radio source with flux of l.lmJy at 20 cm. The spectra of all three quasars show 
strong and broad Lya+NV emission lines, and very strong Lya forest absorption, 
with a mean continuum decrement Da > 0.90. The ARC 3.5m spectrum of SDSSp 
J103027.10+052455.0 shows that over a range of ~ 300A immediately blueward of 
the Lya emission, the average transmitted flux is only 0.003 ± 0.020 times that of the 
continuum level, consistent with zero flux over a ~ 300A range of the Lya forest region, 
and suggesting a tentative detection of the complete Gunn-Peterson trough. The 
existence of strong metal lines in the quasar spectra suggests early metal enrichment 
in the quasar environment. 

The three new objects, together with the previously published ^; = 5.8 quasar 
SDSSp J104433. 04-012502. 2, form a complete color-selected flux-limited sample at 
z ^ 5.8. We estimate the selection function of this sample, taking into account the 
estimated variations in the quasar spectral energy distribution as well as observational 
photometric errors. We find that at ^; = 6, the comoving density of luminous quasars at 
Mum < -26.8 {Hq = 50 km s"^ Mpc^S = 1) is 1.1 x IQ-'^ Mpc"^. This is a factor 
of ~ 2 lower than that at z ^ 5, and is consistent with an extrapolation of the observed 
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quasar evolution at z < 5. Using the current sample, we discuss the constraint on 
the shape of the quasar luminosity function and the implications for the contribution 
of quasars to the ionizing background at z ~ 6. The luminous quasars discussed in 
the paper have central black hole masses of several times lO^M© by the Eddington 
argument, with likely dark halo masses of order IO^^Mq. Their observed space density 
provides a sensitive test of models of quasar and galaxy formation at high redshift. 



Subject headings: quasars:general; quasars: absorption line; quasars: emission line; 
intergalactic medium 



1. Introduction 



High-redshift quasars provide direct probes of the epoch when the first generation of galaxies 
and quasars formed. The absorption spectra of these quasars reveal the state of the intergalactic 
medium (IGM) close to the reionization epoch ( [Haiman Sz Loeb 1999| , Miralda-Escude, Haehnelt 
fc Rees 200d| , [Madau fc Rees 2000| , |Cen fc Haiman 200q) . The lack of a Gunn-Peterson trough 
(Shklovsky 1964, Scheuer 1965, Gunn & Peterson 1965) in the spectrum of the luminous quasar 
SDSSp J1O4433.O4-O125O2.20 at z = 5.8 (SDSS 1044-0125 for brevity, |Fan et al. 20"00^ , for an 
updated redshift of this object, see Poodrich et al. 2001 and Djorgovski et al. 2001 ) indicates 
that the universe was already highly ionized at that redshift. Assuming that SDSS 1044-0125 is 
radiating at the Eddington luminosity, this object contains a central black hole of several billion 
solar masses. The assembly of such massive objects in a timescale shorter than 1 Gyr yields 



constraints on models of the formation of massive black holes (e.g., Haiman fc Loeb 2001 ). The 
abundance and evolution of such quasars can provide sensitive tests for models of quasar and 
galaxy evolution. 



The Sloan Digital Sky Survey (SDSS; York et al. 2000| , Stoughton et al. 2001 ) uses a 
dedicated 2.5m telescope and a large format CCD camera ( punn et al. 1998 ) at the Apache Point 
Observatory in New Mexico to obtain images in five broad bands {u, g, r, i and z ^ centered at 
3551, 4686, 6166, 7480 and 8932 A, respectively; [Fukugita et al. 1996|) over 10,000 deg^ of high 
Galactic latitude sky. Over 200 quasars at z > 3.5 have been discovered to date from the SDSS 
multicolor imaging data (Fan et al. 1999a, 2000a,c, 2001a, [Zheng et al. 200q , Schneider et al. 



2000, 2001, Anderson et al. 2001). The inclusion of the reddest band, z, in principle enables the 
discovery of quasars up to z ~ 6.5 from the SDSS data (Fan et al. 2000c). 



^^The naming convention for the SDSS sources is SDSSp JHHMMSS.SSiDDMMSS.S, where "p" stands for the 
preliminary SDSS astrometry, and the positions are expressed in J2000.0 coordinates. The astrometry is accurate to 
better than 0.2" in each coordinate. 



Following Stoughton et al. 2001, we refer to the SDSS passbands as u, g, r, i and z. As the SDSS photometric 
calibration system is still being finalized, the SDSS photometry presented here is referred to as it*, g* , r* , i* and z* . 
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Encouraged by the discovery of SDSS 1044-0125, we have carried out a systematic survey of 
i-dropout quasars {z ^ 5.8, whose Lyman a forest is entirely in the i band) based on the SDSS 
imaging data. The goals of this survey are: (1) to define a color-selected flux-limited complete 
sample of quasars at z ~ 6 in order to study the evolution of the quasar spatial density at high 
redshift; (2) to provide multiple lines of sight for studying the nature of the high-redshift IGM; (3) 
to study the intrinsic properties of quasars and the metallicity of the quasar environment at high 
redshift. 

In this paper, we present the first results of this survey, covering an area of 1550 deg^. In 
§2, we present our color selection procedures for z > 5.8 quasars. Special software is developed 
to control the number of false single-band detections in the imaging data. Follow-up optical 
and near-infrared photometry is used to separate high-redshift quasar candidates from the more 
numerous cool dwarfs which have similar optical colors. In §3, we describe the discovery of three 
new quasars at z = 5.82, 5.99 and 6.28. We discuss their spectral properties and the implications 
for the metallicity of the quasar environment and the status of the IGM in §4. These three 
quasars, combined with SDSS 1044-0125, form a complete color-selected sample at z* < 20.2. We 
calculate the survey selection function and derive the spatial density of z ~ 6 quasars in §5. In §6, 
we discuss the cosmological implications of the observations. Follow-up spectroscopic observations 
carried out with the Keck telescope will be presented by Becker et al. (2001; Paper II). 

We present our results using two different cosmologies throughout the paper: a A-dominated 
universe with Hq = 65 km s~^ Mpc~^, A = 0.65 and O = 0.35 ( pstriker fc Steinhardt 1995 , Krauss 



fc Turner 1995 ), which is referred to as the A- model; and a Einstein-de Sitter universe with = 1 
and Hq = 50 km s"^ Mpc~^, which we refer to as the ^1 = 1 model in this paper. The age of the 
universe is 13.9 Gyr and 13.0 Gyr in the A-model and the Q = 1 model, respectively. 



2. Candidate Selection 

2.1. Colors of z > 5.8 quasars 

The observed colors of quasars evolve strongly with redshift, as the intrinsic emission and the 
intervening absorption spectral features move through the SDSS filter system ( [Fan 199£ , Richards 



et al. 2001a| ). The strong Lya forest absorption blueward of the Lya emission line enters the SDSS 
g, r and i bands at redshifts of about 3.6, 4.6 and 5.5, respectively. At z ^ 5.7, the Lya emission 
line begins to move out of the SDSS i filter. With an average i* — z* ^ 2, these high-redshift 
quasars have undetectable flux in the bluer u, g and r bands, and little flux (sometimes also 
undetectable) in i band. They become z-dropout objects with only one measurable color in SDSS 
photometry. 



Figure 1 presents the i* — z* vs. z* color-magnitude diagram for 50,000 stellar sources selected 
at random from the SDSS imaging data. It also shows the locations of all the i-dropout objects 
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(including L/T dwarfs and high-redshift quasars, see below) selected in the 1550 deg^ survey 
area. The median track of simulated quasars (sec §5.1) with absolute magnitude M1450 = —27, 
representative of the quasars discussed in this paper, is also plotted. A simple cut at i* — z* > 2.2 
selects quasar candidates at z > 5.8. At 2; > 6.0, the Lyct forest absorption enters the z filter, and 
the z* flux decreases rapidly with redshift with increasing amount of absorption, while i* — z* 
remains approximately constant, as the increasing absorption in the i and z bands are roughly 
equal. At 2; ~ 6.6, the Lya emission line has redshifted to 9240A, well past the effective wavelength 
of the z filter. Most of the z flux is absorbed by the Lya forest, and the object disappears from 
optical images altogether. This is the upper limit to the redshift of objects detectable with SDSS 
imaging data. 

Quasars at z > 5.8 are extremely rare on the sky. We found four i-dropout quasars in 1550 
deg^. For comparison, over the same area, there are about 15 million objects detected in the z 
band above 6-cr, and about 6.5 million cosmic ray hits in the z band. The quasars are faint, with 
low signal-to-noise ratio (S/N) photometry from the SDSS imaging data. Therefore, the key to 
a successful selection process is an efficient elimination of potential contaminants. We face three 
technical challenges: 

(1) Elimination of false z-band only detections, composed primarily of cosmic rays, but also 
including satellite trails, electronic ghost images and bleed trails from bright stars. While the 
recognition of the majority of cosmic rays is relatively straightforward, even a tiny fraction of 
mis-classified cosmic rays in the z band, which would then appear as i-dropout candidates, would 
dominate over the much rarer true quasars. Note that SDSS scans the sky only once in the course 
of the survey except in a small fraction of overlapping area; thus repeat observations cannot be 
used to recognize cosmic rays. Moreover, SDSS uses thick CCDs for the z band. A cosmic ray hit 
on these thick devices typically occupies more than a single pixel, making the separation of real 
detections from cosmic ray hits more difficult, especially when the seeing is good. 

(2) Reliability of z* photometry. As is evident from Figure 1, the i-dropout objects arc 
distributed on the extreme tail of the i* — z* distribution. Under typical seeing, cr(z*) ~ 0.1 for 
objects with z* 20. Objects in the tail of the error distribution could scatter into the selected 
region of color space, and pose a serious contamination when a large number of objects is searched. 

(3) Separating quasars and cool dwarfs. As shown in §3, the surface density of cool dwarfs, 
with spectral types ranging from mid-L to T, which have i* — z* > 2.2 (Geballe et al. 2001, Leggett 
et al. 2001), is ~ 15 times higher than that of z > 5.8 quasars. We use J-band photometry to 
separate these two classes of objects (see below, also Fan et al. 2000b, Zheng et al. 2000). 

2.2. Selection Procedure 

Because of the rarity of high-redshift quasars, and the overwhelming number of contaminants, 
our photometric selection procedure of z > 5.8 quasar candidates is quite involved, and includes 
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five separate steps: 

1. Selection of i-dropout sources from tlie SDSS database. 

2. Running an improved cosmic ray classifier and visually inspecting all remaining candidates 
to reduce the number of false candidates due to cosmic rays. 



3. Matching the candidate list to the Two- Micron All Sky Survey (2MASS, [Skrutskie et 
al. 1997) in the public release area of 2MASS, to obtain JHK photometry. 



4. Independent z photometry to further eliminate false and biased z detections. 

5. J-band photometry of objects not in 2MASS to separate quasar and cool dwarf candidates. 
This then leaves us with a sample for follow-up spectroscopy. 

In this and the subsequent section, we describe this process in detail. The techniques and 
rationale for these various cuts are described in the present section, while the observations 
themselves are described in § |3[ 

2.2.1. Selection of i- dropout objects from SDSS imaging data 
The SDSS photometric data are processed by a series of automated pipelines to carry out 



astrometric and photometric measurements. The photometric pipeline {Photo, Lupton et al. 2001 ) 
reduces the data from the imaging camera and produces corrected images and object catalogs. 
Photo classifies detected objects either as "star" (point source, consistent with the point spread 
function, PSF) or "galaxy" (extended source). However, we select all z-dropout objects irrespective 
of their star-galaxy classification, as the star-galaxy separation becomes less reliable at low S/N. 
Cosmic rays are recognized as objects with profile gradients significantly steeper than that of 
the PSF, and are interpolated over; they are not output explicitly by the pipeline. While Photo 
classifies more than 99.9% of the cosmic rays correctly, visual inspection shows that a majority of 
objects with z-only detections that are classified as "stars" still have the appearance of cosmic 
rays. This is further confirmed by comparing those areas of the sky with multiple observations. 
Therefore, we have recently improved the performance of the object classifier, as follows. 

For every object detected by the imaging pipeline, we count the pixels in an 11 x 11 pixel 
region (4.4" x 4.4") centered on the object, for which 

PSF(d) I{d)+c*N[Iid)] 

PSF(O) ^ I{0)-c*N[I{0)] ^ ' 

where /(O) is the object's central intensity, I{d) is the intensity of the pixel in question, PSF(O) 
and PSF(d) are the equivalent quantities for the locally-determined model of the PSF, the noise 
is the photon noise with an additional 0.05/(0) added in quadrature, and c = 1.5. This criterion 
searches for objects with profile gradients significantly sharper than that of the PSF. If more than 
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3 pixels satisfy this condition, we consider the object to be a candidate cosmic ray, and reject it. 

The photometric pipehne assigns detailed flags for each detected object, indicating objects 
whose photometry (and therefore, colors) may be in error. We reject objects flagged as being 
saturated in any band, lying on the bleed trail of a saturated star, or overlapping the edge of 
the image boundary. We pay close attention to objects blended with close neighbors, and reject 
deblended children with PEAKCENTER, NOTCHECKED, or DEBLENDJJOPEAK set in the i and z bands, 
which are typical signs of deblending problems. For a complete description of the definition of 
flags, see ^toughton et al. 2001 . The flag-checking procedure used here is similar to that used in 



the SDSS quasar target selection pipeline ( Richards et al. 2001b ) 



We have re-classified all the objects with i* — z* > 2.2 using the new cosmic ray classifier. We 
found that depending on the seeing of a specific run, 50% to 70% of the z-dropout objects formerly 
classified as "stars" are now classified as cosmic rays, while a negligible fraction of the "real" 
objects (based on multiple observations and objects with confirmed spectra) are mis-classified as 
"cosmic rays". We then visually inspect the z images of the remaining candidates, and reject a 
further 25% as cosmic rays and other artifacts. 

After the improved cosmic ray rejection and visual inspection, we carry out independent 
follow-up z photometry of all the remaining candidates (of which there are 121; see below) using 
the Apache Point 3.5m and other telescopes (see §3.1). This serves two purposes: (a) to confirm 
that the objects are indeed real; and (b) to measure the z* magnitude more accurately. Since 
our z-dropout candidates have the reddest i* — z* colors, we tend to pick up objects lying on the 
several a tail of the z* error distribution (note that the z-band only objects are very faint, and 
therefore are measured at low S/N); they will thus have artificially red i* — z* color and blue 
z* — J color, and are more likely to be selected as z > 5.8 quasar candidates (see below). To 
include those objects in the sample would lower our spectroscopic success rate considerably. 



2.2.2. Near-IR Photometry: separating quasars and L/T dwarfs 

As discussed by Fan et al (2000b) and Zheng et al. (2000), and illustrated in Figure 1, the 
major astrophysical contaminants of i-dropout quasars are very cool dwarfs with spectral types 
L and T (Strauss et al. 1999, Tsvetanov et al. 2000, Fan et al. 2000a, Leggett et al. 2000, 2001, 
Geballe et al. 2001). These have effective temperatures T'e// < ISOOX, and i* — z* > 2, and cannot 
be distinguished from z > 5.8 quasars using SDSS optical photometry alone. Figure 2 presents the 
i* — z* vs. z* — J color-color diagram of the i-dropout sample. Finlator et al. (2000) have matched 
the SDSS photometric catalog in ~ 50 deg^ with that of 2MASS in the JHK bands. The small 
dots in Figure 2 show the SDSS-2MASS matches and indicate the location of normal stars on this 
color-color diagram. The track of quasar colors shows that by using z* — J < 1.5, i-dropout cool 
dwarfs and high-redshift quasars are clearly separated, as the z* — J color of high-redshift quasars 
is dominated by the blue power law continuum, while the cool temperature of L/T dwarfs results 
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in very red z* — J color, typically larger than 2. Therefore we determine whether the candidate has 
been detected by 2MASS. While the non-detection of a faint candidate by 2MASS does not give a 
tight constraint on the spectral shape of the object, most of the brighter objects {z* < 19) that 
are located in the area of the 2MASS publicly released data are indeed detected by 2MASS, with 
z* — J ^ 2, indicating that most of them are cool dwarfs. For those not matched with 2MASS, we 
carry out follow-up J-band photometry (see §3.1). 



2.2.3. Color selection criteria 
The final photometric selection criteria for this survey are: 

(a) z* < 20.2, 
(6) a{z) < 0.1, 

(c) i* -z* > 2.2, ^ ^ 

(d) z* - J < 1.5. 

We also require that the object not be detected in any of the other three bands u, g and r. 

The magnitudes and colors in Eq. (1) are not corrected for interstellar extinction. This is a 
small effect, which has been taken into account in the selection function calculation (§5.1). In the 
selection, only objects with S/N > 10 in the z band are considered (item [b]), as the i — z color 
becomes very unreliable for the faintest sources and would have selected many more candidates 
with a more inclusive cut. Item (c) selects quasars with z ^ 5.8, and item (d) separates quasars 
from cool dwarfs. Note that the SDSS photometry is on the AB magnitude system (Fukugita 
et al. 1995), while the J-band photometry is on the Vega-based system. The selected area is 
illustrated as the shaded region in Figure 2. 

3. Discovery of Three New Quasars 

3.1. Photometric Observations 

Table 1 lists the details of the 25 SDSS photometric runs used in this survey. The projection 
of these runs on the sky is illustrated in Figure 3. The entire region is located in the Northern 
Galactic Cap. In Table 1, the location of each run is presented in terms of survey stripe number 
(n) and the range of survey great circle longitude (Amin and Amax)- When taking data, the SDSS 
telescope moves along a series of great circles on the sky and the photometric camera drift-scans 
at the sidereal rate, with an eff'ective exposure time of 54.1 seconds. The survey coordinate system 
(A, 7]) is a spherical system with poles at a2ooo = 95°, and 275°, ^2000 = 0°. A = 0°, ?? = 0° is 
located at 02000 = 185°, ^2000 = 35° (York et al. 2000, [Stoughton et al. 2001| ). Each drift scan 
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tracks a survey stripe n centered on a constant 77, given by 

?7 = (n - 10) X 2.5° - 32.5. (3) 

Two scans, or strips, one offset to the north and one to the south, are required to fill a stripe. 
The 25 SDSS imaging runs used in the i-dropout survey in this paper were taken between March 
20, 1999 (run 745) and March 19, 2001 (run 2190) over a range of seeing conditions. The seeing 
of SDSS images is characterized by the psf Width parameter, defined as the effective width of the 
best fit double-Gaussian PSF model. For a single Gaussian profile, psf Width = 1.06 FWHM. 
Only those runs, or parts of runs, with psf Width in z band in the 4th column of the SDSS camera 
(near the middle of the focal plane) better than 1.8" are used in this paper. The median seeing 
values (over all six columns of the camera) in the i and z bands arc also listed in Table 1. Due to 
the complicated geometry of the survey and overlap between runs, we estimated the total survey 
area by Monte-Carlo integration - calculating the fraction of randomly distributed points on the 
celestial sphere falling into the boundary of any of the runs listed in Table 1. The total area 
covered is 1550 dcg^, excluding regions with poor seeing. The median psf Width is ^ 1.5" in the i 
and z bands over the whole area. 

Following the procedures described in §2, we first search the SDSS database and eliminate 
false detections (Steps 1 and 2) based on the SDSS images. This results in a total of 121 objects 
with i* — z* > 2.2 that require further z and J band photometry. The candidate list includes five 
known T dwarfs from the SDSS (Strauss et al. 1999, Tsvetanov et al. 2000, Leggett et al. 2000), 
one known T dwarf from 2MASS (Burgasser et al. 1999), and the known z = 5.8 quasar SDSS 
1044-0125. Thirty-five objects are located in the 2MASS first and second incremental data release 
area. Of these, 13 objects are detected in 2MASS JHK photometry and have colors of L to T 
dwarfs. 

Independent z photometry (Step 4 in §2.2) was carried out using the Seaver Prototype 
Imaging camera (SPICAM) in the SDSS z filter on the ARC 3.5m telescope at the Apache 
Point Observatory in several nights between February and April 2001. SPICAM has a backside 
illuminated SITe 2048 x 2048 pixel CCD with a field of view of 4.8'. The pixel scale is 0.14", 
and the typical seeing was 1". The quantum efficiency peaks at 89% at 6500A, and is about 
49% at 9OOOA. The exposure times were between 90 and 180 seconds, depending on the weather 
conditions. Note that since the SDSS photometry itself provides many local reference stars, 
photometric weather is not required for z photometry. The S/N is typically twice that of the 
SDSS imaging. The SPICAM images arc reduced using normal IRAF procedures. Sky flats are 
constructed from the scaled median of target images, z* magnitudes are measured with aperture 
photometry with a radius of ~ 2", and calibrated by local SDSS standards. 

Follow-up J band infrared photometry (Step 5 in §2.2) was carried out on photometric nights 
using several telescopes, between February and May 2001. The majority of the targets were 
observed with the GRIM II instrument (the near infrared GRIsm spectrometer and IMager), also 
on the ARC 3.5m. It uses a 256 x 256 HgCdTe detector and covers a field of view of 1' at f/10. 
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JHK standards from Persson et al. (1998) were monitored throughout the night. We used similar 
IRAF procedures to reduce the GRIM II data. A small number of objects were observed with 
the UFTI instrument, a near IR camera with a 1024 x 1024 HgCdTc detector covering 92" x 92" 
on the 3.8m United Kingdom Infrared Telescope (UKIRT) on Mauna Kca, and with the MAGIC 
instrument, a near-infrared imager and low-resolution infrared spectrograph at the 2.2m telescope 
at Calar Alto Observatory in Spain. The typical photometric error of our infrared photometry 
is 0.05 - 0.10 mag, and is dominated by calibration errors. Since quasars and L/T dwarfs have 
z* ~ J colors that differ by more than 0.5 mag, this error has very little effect on the selection 
efficiency (§5.2). 

We attempted to follow the sequence described in §2.2 to select z-dropout quasar candidates. 
But as we had to use photometric nights for the IR photometry, we changed the observing sequence 
in two cases: (a) on a dry photometric night, IR photometry had the highest priority; (b) when 
the weather was non-photometric, for objects that had SPICAM z photometry but did not yet 
have J band photometry, we carried out spectroscopy directly, without pre-selection using z* — J 
color. In the latter case, we started from either the brightest or the reddest object, depending on 
weather. It is worth noting that the rapid instrument change mechanism at the ARC 3.5m 
15 minutes) and the excellent weather monitoring system at APO allows a very efficient use of 
telescope time. In the end, we obtained photometry and/or spectroscopy of our full sample. 

3.2. Spectroscopic Observations 

Table 2 summarizes the classifications of the i-dropout sample. Using 2MASS matches 
and follow-up z and J photometry, we found that among the 121 z-dropout candidates, 35 are 
false detections, most likely cosmic rays that still passed the first two steps of selection; 71 have 
z* — J > 2, and are classified as cool dwarfs. Ten objects are classified as cool dwarfs based on 
spectroscopy; so no J photometry was necessary. The spectra of these dwarfs were obtained using 
the Double Imaging Spectrograph (DIS) on the ARC 3.5m between March and May 2001, and 
using the Echelle Spectrograph and Imager (ESI) on Keck II telescope in March 2001. One of 
the candidates, SDSSp J081948.97-F420930.2, has an intermediate z* - J color {z* - J = 1.7); it 
is a low-ionization BAL quasar at z = 2.05 (see Becker et al. 1997, Hall et al. 2001). The Keck 
spectrum of this object is presented in Hall et al. (2001). Finally, four objects have i* — z* > 2.2 
and z* — J < 1.5, including the previously known z = 5.8 quasar SDSS 1044-0125, and three new 
quasars at z > 5.8. Note that in this survey, all objects that satisfy the color selection criteria 
(Eq.[l]) indeed turn out to be z ;^ 5.8 quasars from their spectra. The magnitude and colors of 
these objects are plotted on Figures 1 and 2. 

Among the cool dwarfs found in this survey, eleven are T dwarfs. They form a complete 
fiux-limited T dwarf sample. Six of these objects were previously known; the optical-IR spectra of 
four more of them are presented in Geballe et al. (2001). In a separate paper, we will present the 
photometry of L and T dwarfs in the z-dropout sample, the observations of the final T dwarf, and 
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the analysis of the surface density and spatial densities of T dwarfs in the solar neighborhood. 

The optical spectra of the three new quasars were obtained using the DIS on the ARC 3.5m 
between March and May 2001. The instrument and data reduction procedure are described 
in detail by Fan et al. (1999a). Figure 4 presents the z band finding charts of the three 
new quasars, SDSSp J083643.85+005453.3 {z = 5.82, SDSS 0836+0054 for brevity), SDSSp 
J130608.26+035636.3 (SDSS 1306+0356, z = 5.99) and SDSSp J103027.10+052455.0 (SDSS 
1030+0524, z = 6.28). Table 3 presents their photometric properties, and Figure 5 shows the ARC 
3.5m discovery spectra of the three objects. The exposure time of each spectrum is 3600 seconds. 
The spectra in Figure 5 have been smoothed to a resolution of ~ 20A. The flux calibration has 
been adjusted to match the SDSS z* magnitudes. The locations of the Lya and NV emission lines 
are indicated in the Figure. The expected positions of Ly/3+OVI and the Lyman Limit are also 
shown. 

A near-IR spectrum of SDSS 1030+0524 was obtained on the night of 11 June 2001, 
using the NIRSPEC ( [McLean et al. 1998 ) instrument on the Keck II telescope, as part of the 



Gemini/NIRSPEC service observing program, carried out by Gemini staff scientists Tom Geballe 
and Marianne Takamiya. Observations were taken in grating mode, with the N2 blocking filter 
and a 0.76" slit, which gives a wavelength coverage of 1.09 /xm - 1.29 fim, and a spectral resolution 
of ~ 1500. A total of 10 exposures, each with 300 sec exposure time, were taken using an ABBA 
sequence, nodding along the slit. The weather was photometric with ~ 0.7" seeing. HIP 48414, 
a bright AlV star, was observed as an atmospheric calibration star at similar airmass to that of 
the quasar observation. The NIRSPEC data were reduced using the WMKONSPEC package, an 
IRAF package developed at Keck Observatory to reduce NIRSPEC data. The flux calibration has 
been adjusted to match its measured J-band magnitude. The final reduced spectrum, binned to 
4Apixel~^, is shown in Figure 6. For comparison, the optical spectrum of SDSS 1030+0524 is also 
shown, and the locations of the Lya, NV and CIV emission lines are indicated. 



4. Spectral Properties 

The discovery spectra of the three objects show unambiguous signatures of very high-redshift 
quasars, similar to that of SDSS 1044-0125: strong, broad and asymmetric Lya+NV emission 
lines, with sharp discontinuities to the blue side, due to the onset of very strong Lya absorption. 
The discovery spectra are of relatively low S/N. Accurate redshift determination from these 
spectra are difficult. The Lya emission line is severely affected by the Lya forest absorption, and 
the NV line is blended with Lya. The low S/N does not allow us to use weaker lines such as 
OI+SillA 1302 or SiIV+OIV]A 1400 to determine an accurate redshift. The next strong emission 
line in the quasar spectrum is CIVA 1549, which is located beyond 1 /im for z > 5.5, and requires 
near-IR spectroscopy. Therefore, we use the near IR spectrum to determine the redshift of SDSS 
1030+0524. The redshifts of the other two quasars were determined using high S/N Keck spectra 
described in Paper II. The uncertainties in all these redshifts are of order 0.02. In this section, we 
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first briefly comment on the major spectral signatures of each quasar (§4.1). We use the emission 
lines to constrain the metallicity of the quasar environment (§4.2). Then we calculate the average 
absorption in the Lya forest region, based on the discovery spectra, and comment on the possible 
detection of the Gunn-Peterson trough in the spectrum of SDSS 1030+0524 (§4.3). 



4.1. Notes on Individual Objects 

SDSSp J083643.85+005453.3 (z = 5.82). This quasar has a similar redshift to SDSS 
1044-0125. The Lya emission line is very broad and strong. A separate NVA1240 component is 
tentatively detected. The rest frame equivalent width (EW) of the Lya+NV line is ~ 70 A, quite 
typical of lower-redshift quasars (Fan et al. 2001b). The Ly/?+OVI emission line is also clearly 
detected at ~ 7000A; this line was also seen in the spectrum of SDSS 1044-0125. It is evident that 
there is detectable flux blueward of the Lya emission line. The universe is still highly ionized at 
z ~ 5.8. 

SDSS 0836+0054 is an extremely luminous object. With z* = 18.74 and J = 17.89, it has an 
absolute magnitude M1450 = —27.62 (Q = 1 model, see §5.1), and is the most luminous quasar 
discovered at 2; > 4.5 to date; in particular, it is 0.5 mag more luminous than SDSS 1044-1025. 



SDSS 0836+0054 has a radio counterpart in the FIRST radio survey (Becker, White & Helfand 



1995) at 20cm. The radio source is unresolved with a total flux of 1.11 ± 0.15 mJy, and the 
positional match is better than 1". This is the highest-redshift radio loud quasar known, and the 
second radio-loud quasar detected at 2: > 5 (the other one is SDSSp J091316. 56+591921. 5, at 
z = 5.11, with a 20cm flux of 18.1 mJy; see Anderson et al. 2001). 

SDSSp J130608. 26+035636. 3 {z = 5.99). This is the second highest-redshift quasar known 
to date. The Lya emission line is centered at ~ 8500A, and a separate NV component is also 
detected. The rest frame equivalent width of Lya+NV is ~ 60 A. The Lya forest is stronger than 
that in SDSS 0836+0054, but there is still detectable flux in the Lya forest region, immediately 
blueward of the Lya emission line. The apparent absorption at ~ 9400A is not real, but is due to 
imperfect telluric absorption correction in this noisy spectrum. The object has z* = 19.47 and is 
detected at the 4-0" level in the SDSS i band. It is also a luminous object, with M1450 = —26.93 in 
the = 1 model. 

SDSSp J103027.10+052455.0 {z = 6.28). This is the highest-redshift quasar known to 
date. At z* = 20.05 and J = 18.87, it is also the faintest quasar found in the current survey, and 
has M1450 = —26.89 in the Q = 1 model. It is close to the flux limit of the survey, and is only 
detected at lO-cr in z. The object is very red, and is completely undetected in the SDSS i band. 
The discovery spectrum shows very strong Lya+NV emission, with a rest-frame equivalent width 
of ~ 70A, and a tentative detection of a separate NV component. The Lya absorption is strikingly 
strong, with almost no detectable flux blueward of the Lya emission. A strong CIV emission 
line is detected in the near IR spectrum (Figure 6). A Gaussian fit to the CIV line profile gives 
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Zeni = 6.28 ± 0.02, with a rest frame EW = 31.5 ± 8.6A. The continuum shape at 9000 A < A < 
13000 A is consistent with f,y oc z/~'^'^, as indicated by the dashed hue in Figure 6. 



4.2. Metallicity of z ~ 6 Quasars 

Emission Hne ratios can be used to measure the metalhcity of the gas in the Broad Emission 
Line Region (BELR). There is growing evidence from those measurements that BELR have 
roughly solar or higher metallicities even out to z > 4 (e.g., Hamann Sz Ferland 1999| ). We use the 



spectra of our quasars to estimate the quasar metallicity at z ~ 6. 

Of various line ratios, the NV A1240/CIV A1549 ratio and NV A1240/HeII A1640 ratio 
are particularly useful abundance diagnostics ( [Hamann fc Ferland 1993| , Pamann fc Ferland 



19991) . Hamann & Ferland (1993) examined the chemical evolution of BELR gas by applying 
spectral synthesis and chemical enrichment models to the NV/HeH and NV/CIV line ratios. 
They found that nominal BELR parameters predict NV/HeH near unity for solar abundance, 
increasing to ~ 10 for Z ~ WZq. The sensitivity of NV/CIV to metallicity is due to the fact 
that N is a secondary element, with abundance proportional to Z^, while C is a primary element 
with abundance proportional to Z. Their calculations show that for the ISM associated with a 
population with solar abundance, NV/CIV < 0.1. Luminous quasars at 2 < z < 5 show NV/CIV 
~ 0.1 — 2, indicating a metallicity of Zq ^ Z ^ IOZq. 

From the Keck/NIRSPEC J-band spectrum of SDSS 1030+0524, we find that the rest frame 
EW of CIV is 31.5 ± 8.6 A. We do not detect Hell in emission at ~ 11900A, which indicates 
that the Hell EW is smaller than ~ 5A at 3-o". The discovery spectra in Figure 5 show tentative 
detections of NV emission in all three quasars. The quality of the spectra does not allow an 
accurate fit to the Lya+NV line profile in order to derive the equivalent width of NV. If we assume 
the EW of NV to be of the order lOA, and a power law continuum with f,y oc u~^'^, we find the 
flux ratios NV/CIV ^ 0.4, and NV/Hell ^ 3.0. Both limits would imply super-solar metallicity of 
the BELR region. In the calculation of Hamann & Ferland (1993), these line ratios are consistent 
with the "Giant Elliptical" model with fast stellar evolution and a top-heavy initial mass function, 
with 3Zq < Z < WZq. A high resolution, high S/N spectrum is needed to fit the Lya+NV line 
profile and put better constraints on the line ratios (Paper II). However, note that Krolik &; Voit 
(1998) point out that NV can be excited by resonance scattering from the red wing of the Lya 
emission line, thus biasing these metallicity calculations. 



4.3. Average absorption in the Lya forest region 

The luminous high-redshift quasars discussed in this paper are ideal targets for detailed 
studies of the intergalactic medium at high redshift. We calculate average absorption in the Lya 
forest region based on the discovery spectra. The results are summarized in Table 4, which also 
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includes the measurements for SDSS 1044-0125. Three quantities are calculated for each quasar. 
We first estimate the average continuum decrements as: Da,b = ~ f°^^/fv°^^7 where f°^'^ 
and f^°^ are the observed and the unabsorbed continuum fluxes of the quasar, and and Db 
measure the decrements in the region between rest-frame Lya and Ly/? (A = 1050 — 1170 A) and 
between Ly/? and the Lyman Limit (A = 920 — 1050 A), respectively ( Qke fc Korycansky 19821 ). 



Following Fan et al. (2000b), we measure Da and Db assuming a power law continuum z/" with 
a = —0.5. The errors on Da and Db are dominated by the continuum determination. The error 
bars in Table 4 reflect the range of Da and Db-, allowing a to range from —1.5 to +0.5. 

The range of IGM redshifts covered by the window over which Da is calculated is rather large 
(5.0 < Zabs < 5.7 for a quasar at z = 6.0). Therefore we also calculate the transmitted flux ratio 
over a smaller window (of 0.2) centered on a redshift Zabs = -^em — 0.3: 

n^abs) = (f^'yfr) , (l + ^abs - O.l) X 1216A < A < (1 + Zabs + 0.1) X 1216A. (4) 



We assume a power law continuum of f,y oc The error bars in Table 4 reflect only the photon 

noise of the spectra. The quantity T is more noisy than Da, but can be calculated up to within 
a redshift of 0.2 of the quasar; closer than this, it would be affected both by the proximity effect 
from the quasar itself and the blue wing of the Lya emission line (Az ~ 0.2, Oke & Korycansky 
1982). In Table 4 and Figure 7, we show the T(zabs) values for the redshifts Zabs ~ -^em — 0.3. 

From Table 4 and Figure 7, it is evident that the amount of Lya absorption increases rapidly 
with redshift, with Da > 0.9 at Zem > 5.8. For comparison, Da ~ 0.75 at Zem ~ 5 (Songaila et 
al. 1999). The most striking feature, however, is that based on the spectrum of SDSS 1030+0524, 
T{zahs = 6.0) = 0.003 ± 0.020, consistent with no flux detected at all. The flux level is consistent 
with zero from 8400A to 8700A in the spectrum (Figure 5), corresponding to 5.9 < Zabs < 6.15. 
The 1-0" lower limit on the flux decrement in this quasar is a factor of 50. At slightly lower 
redshift, we find f / fconizabs = 5.5 — 5.7) = 0.07 — 0.10 from the other three z > 5.8 quasars in 
the sample. The flux decrement observed in SDSS 1030+0524 is at least a factor of 4 larger than 
that seen in any other quasar at z > 5. There seems to be a drastic change in the amount of Lya 
absorption at z 6. 



4.4. Gunn-Peterson Effect in the Spectrum of SDSS 1030+0524 

The absence of flux over a 300A region in the ARC 3.5m discovery spectrum indicates a 
possible first detection of the complete Gunn-Peterson trough. This detection is still highly 
tentative, due to the uncertainty in sky subtraction in these low S/N spectra. Note that spectra 
taken on different nights show the same non-detection of flux in the Lya forest region. This result 
shows that the fraction of neutral hydrogen has increased substantially between z = 5.7 and z = 6, 
and, if confirmed by high S/N spectroscopy, would show that the universe is approaching the 
epoch of reionization. In Paper II, we use the high-resolution Keck spectrum of SDSS 1030+0524 
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to place stronger constraints on the status of the IGM at z ~ 6, and discuss its imphcation on the 
reionization of the universe in more detail. 

For a uniformly distributed IGM, the Gunn-Peterson (1965) optical depth is: 

If the IGM were mostly neutral, ~ 1, and the Gunn-Peterson opacity would be ~ 10^. 
Even a tiny fraction of neutral hydrogen (of the order of 10~^) in the IGM could result in a 
large optical depth, and undetectable flux in the Lya forest region. Therefore, the existence of 
the Gunn-Peterson trough by itself does not indicate that the object is observed prior to the 
reionization epoch, as it takes only a tiny fractional abundance of neutral hydrogen to create a 
trough. Indeed, both semi-analytic models (Miralda-Escude et al. 2000) and hydro-dynamical 
simulations ( Gnedin 2000D of structure formation show that the Gunn-Peterson trough should 



begin to appear at z ^ 6, while the epoch of overlap of HII regions could have occurred at a higher 
redshift. 



5. Spatial Density of z ~ 6 Quasars 

5.1. A Color-selected Complete Sample at z > 5.8 

The three new z > 5.8 quasars presented in this paper, plus SDSS 1044-0125, comprise a 
complete color-selected flux-limited sample of z > 5.8 quasars at z* < 20.2. The color selection 
criteria are listed in Eq. (1), and the total sky coverage is 1550 deg^, using the SDSS photometric 
runs in Table 1. The continuum properties of this sample are given in Table 5. Following 
Fan et al. (2000c), the quantity ABi2so is defined as the AB magnitude of the continuum at 
rest-frame 1280A, after correcting for interstellar extinction using the map of Schlegel, Finkbeiner 



fe Davis (1998). We extrapolate the continuum to rest-frame 1450A, assuming a continuum shape 



fi, oc i>~^-^ to calculate ABi^^q. Table 5 also lists the absolute magnitudes M1450 and M1280 in 
both the A and Q = 1 models. In the next subsections, we first calculate the selection function of 
our color-selected sample, then derive the total spatial density of luminous quasars at z ~ 6. 



5.2. Selection Function 

The selection function p(Mi45o, z) is defined as the probability that a quasar of a given M1450 
and z will satisfy the selection criteria (eq. [1]). We calculate this using a Monte-Carlo simulation 
of quasar colors, based on the quasar spectral model described in Fan (1999). We follow the 
procedures described by Fan et al. (2001a), with several modifications: 

(1) J band magnitudes are included in the calculation. We use the J band filter curve 
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measured for the filters used in the GRIM II instrument on the ARC 3.5m (we assume that the 
the UFTI and MAGIC instruments have similar J band filter curves). It is very similar to the 
filter used in the 2MASS survey. The detailed filter curve is available from the authors. Since 
the J band flux of the quasar is dominated by the power law continuum, any small differences 
between the J band filter curves of different instruments has little effect on the selection function. 

(2) SDSS photometric errors, as a function of seeing, are modeled in more detail. The 
selection criteria include (t{z*) < 0.10, which is a strong function of the z* magnitude and the 
observed seeing. We first construct the distribution of i and z seeing for the entire survey area. 
In the simulation, we randomly draw the i and z seeings from this distribution and use them to 
calculate the photometric error terms of the SDSS magnitudes. 

(3) The selection is based on the z* magnitudes without correcting for interstellar extinction. 
We therefore construct the E{B — V) distribution of the entire survey area based on the dust map 
of Schlegel, Finkbeiner fc Davis (1998 J , and assign values of E{B — V) drawn at random from this 



distribution, to the quasars in the simulation. 

Following Fan et al. (2001a), we model the quasar intrinsic spectral energy distributions 
(SEDs) using a power law continuum f^, oc with a = 0.79 ± 0.34, plus a series of broad 
emission lines, with the rest-frame equivalent width of Lya+NV with a mean of 69.3A and a 
standard deviation of 18. oA. We use the quasar absorption line models in Fan (1999) to simulate 
the Lya absorption, which is critical to predicting the colors of high-redshift quasars. In this 
model, we assume that the number density of Lya forest lines evolves as N{z) oc (1 + z)^-^. 

The results of the selection function are summarized in Figure 8 for both cosmologies. 
The heavy lines in the figure represent the 5% contour, indicating the survey limit in absolute 
magnitude at each redshift. From the figure, we notice: 

(1) The survey is sensitive to the redshift range 5.8 ^ z ^ 6.3. The lower limit is due to the 
cut i* — z* > 2.2. The upper limit is due to both the cut z* — J > 1.5, and the fact that quasars 
begin to disappear from z band images because of Lya absorption in the z band (§2.1). 

(2) At 5.8 ^ z ^ 6.0, the survey reaches M1450 ~ —26.5, and is rather complete {p > 0.8) 
at the bright end. At z > 6.0, it becomes increasingly incomplete even for M1450 < —27. The 
survey limiting absolute magnitude as a function of redshift changes much faster than the change 
in distance modulus would indicate, due to increasing Lya absorption in the z band. 

Figure 8 also plots the locations of the four quasars in the sample. Three of them are located 
in the region of high selection probability (p > 0.6), while the highest-redshift and faintest object, 
SDSS 1030+0524, has a selection probability of only ~ 20%. 

These calculations are based on the statistics of the quasar emission line, continuum and 
absorption systems observed at much lower redshift, z ~ 4. To test the sensitivity of our 
conclusions to our imperfect knowledge of quasar SEDs at z ~ 6, we repeated the calculations 
with various continuum slopes and emission line strengths. If we assume a different slope for 
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the power-law continuum, (a) = 0.3 rather than 0.79, we find that the selection probabilities 
of the quasars change by an average of 3%; similarly, when we change the average equivalent 
width of the Lya+NV blend from 69A to 39A, the selection probabilities change by 2%. Thus, 
given the small size of the quasar sample, the error in the estimated spatial density (§5.3 below) 
is dominated by small number statistics rather than the assumptions made in the quasar color 
simulation. However, using these simulations we find that the selection criteria in Eq. (1) will 
not select quasars with undetectable emission lines (EW(Lya) ~ 0, e.g., SDSS 1533-0039, [Fan et| 



al. 1999b ). Such objects will have bluer i — z and redder z — J colors, and are located much closer 
to the stellar locus. They can, however, be selected using a more relaxed color cut. e.g., using 
i — z > 2.0 and z — J < 1.8, the average selection probability of these unusual quasars will be 
higher than 40% for bright quasars at z < 6.2. 

The i* — z* color of z ~ 6 quasars is mainly determined by the Lya absorption. We have very 
little idea of how the Lya forest evolves with redshift at z 6. However, as shown in §4.3, there is 
tentative evidence that the evolution of average absorption is stronger than if we assume a simple 
extrapolation from lower redshift. If that is the case, then at a given redshift, we will overestimate 
the i flux in our models, while the z flux remains almost unchanged, as it is dominated by the 
quasar intrinsic spectrum (continuum -|- emission lines) redward of Lya. Thus the true i* — z* 
color of quasars is redder than the model we assume here, while the z* — J color changes little. 
This means that our i* — z* color cut would be sensitive to quasars at redshifts slightly below 5.8; 
the selection probability at z > 5.8 would remain unchanged. 

The presence of foreground dust in the quasar environment would redden the quasar colors. 
In Figure 2, we indicate the reddening vector for E{B — V) = 0.10 for a quasar at z = 6. As 
discussed in Fan et al. (2001b), the number density of high- redshift quasars could be significantly 
affected by the uncertain amount of dust extinction. 



5.3. Spatial Density of Luminous Quasars at z ~ 6 



With four quasars in the sample, we can only derive the total spatial density of the quasars 
in the redshift and luminosity range the survey covers. We calculate this quantity using the 1/Va 
method, following the discussion in Fan et al. (2001a). For each quasar, the volume over which it 
would have been observed in our survey is 



Va= [ p{Mi45o,z)^dz, 
Jaz dz 



(6) 



where the integral extends over the redshift range 5.7 < z < 6.6. The total spatial density and its 
statistical uncertainty can be estimated as: 
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Using the selection function presented in §5.1, we find that at the average redshift 
of {z) = 5.97, p{M < -26.8) = (1.14 ± 0.58) x 10"^ Mpc'^^ in the O = 1 model, and 
p{M < -27.1) = (0.70 ± 0.35) x 10"^ Mpc'^ in the A-model. 

Schmidt, Schneider & Gunn (1995, hereafter SSG) derive the high-rcdshift quasar luminosity 
function in the range 2.7 < z < 4.75 using 90 quasars at M1450 <J —26. Fan et al. (2001b) 
calculate the evolution of the quasar luminosity function over the range 3.6 < z < 5.0 and 
—27.5 < M1450 < —25.5, using a sample of 39 quasars; the two luminosity functions are in good 
agreement. In Figure 9, we show the density of quasars at M1450 < —26.8 found in this paper, 
along with the results from SSG, Fan et al. (2001b), and the 2dF survey (Boyle et al. 2000) at 
z < 2.5. The quasar density at z ~ 6 found in this paper is consistent with extrapolating the 
best-fit quasar luminosity functions in the range 3 ^ z ^ 5 from SSG and Fan et al. (2001b). 
Both SSG and Fan et al. (2001b) assume a single power law luminosity function at the bright 
end and that the number density of quasars declines exponentially as a function of redshift. For 
example, in the Q = 1 model, the maximum likelihood results of Fan et al. (2001b) find that at 
z = 5, p(Mi45o < -26.8) = 1.98 x 10~^ Mpc"^. Extrapolating the result to z = 6, we predict 
/9(Mi45o < —26.8) = 0.65 x 10~^ Mpc~^, within 1-a of the spatial density estimated in this paper. 
From z = 5 to z = 6, the density of luminous quasars drops by a factor of ^ 1.8. This drop is 
consistent with a decline of a factor of 3 per unit redshift found in Fan et al. (2001b) and 2.7 per 
unit redshift found in SSG. 



6. Discussion 

At M1450 —27, the luminous quasars described in this paper are likely to reside in the 
most massive systems at high redshift. We estimate the masses of the central black holes in these 
quasars following the assumptions of Fan et al. (2000c) : the quasars arc emitting at the Eddington 
luminosity with SED modeled by Elvis et al. (1994), and the fluxes of the observed quasars are 
not significantly magnified by gravitational lensing or beaming. We find that in the A-modcl, the 
estimated black hole masses are Mbh = 4.8 x lO^M©, 2.0 x lO^M©, and 1.9 x lO^M© for SDSS 
0836+0054, SDSS 1306+0356 and SDSS 1030+0524, respectively Under the same assumptions, 
Mbh = 3.4 x IO^Mq for SDSS 1044-0125. In the = 1 model, the estimated black hole masses 
are about 20% lower. Note that the Elvis ct al. (1994) quasar SED, which we use to estimate 
the (substantial) bolometric correction from the flux at rest wavelength of 1450A, is based on 
quasars at z <^ 2. The bolometric correction of quasars at z 6 could in fact be quite different, 
but to determine this will require observations of these high-redshift quasars over a large range of 
wavelengths. 

There is no established relation between black hole mass and the mass of the galactic 
bulge in which it resides, at high redshift. If we make the rather large assumption that the 
Magorrian et al. (1998) relation, MBH/Mi^ige = 3 x 10~^, determined at 2; ~ 0, also holds at high 
redshift, we find that in the A-model, M^uige = (6 — 16) x lO^-'^M© for the quasars in this paper. 
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Assuming further that the ratio of the bulge to dark matter halo mass is of order the ratio of ^If, 
to 0,M, we find that these quasars reside in dark matter halos with mass M^aio ~ lO^^M© for 
Ofj/r^M ^ 10. However, Kauffmann & Haehnelt (2000) argue that the ratio of black hole mass to 
bulge mass could be smaller at high-redshift (see also |Ridgway et al. 2001 ). Alternatively, using 



the Gebhardt et al. (2000) z ~ relation between black hole mass and bulge velocity dispersion 
Mbh = 1-2 X lO*^M0(cJbuige/2OO km s~''f-'^^, we find cJbuige = (420 - 530) kms"^ 

Note that we assume that these four high redshift quasars are not magnified by gravitational 
lensing. In a = 1 universe, the lensing probability rises steeply with redshift up to z ~ 1, 
then becomes rather flat. In a A-dominated model, the increase in lensing probability stays 
steep to somewhat higher redshift, but flattens above the redshift at which the universe is 
matter-dominated. Using the equations of Turner (1990), we find the raw lensing probability to 
be ~ 0.01 at z ~ 6 for reasonable assumptions for the lensing optical depth. The probability of 
magnification bias, however, could be much larger since the sample in this paper is a flux-limited 
sample, and the observed quasars are on the tail of a steep luminosity distribution (Blandford &; 
Narayan 1992). A detailed calculation of magnification bias is beyond the scope of this paper. 

Following Turner (1991) and Haiman & Loeb (2001), we consider the growth of supermassive 
black holes in high-redshift quasars. The e-folding timescale for black hole growth is 4 x 10^(q^)77^^ 
yr, where t/~^ is the ratio of the quasar bolometric luminosity to the Eddington luminosity, and 
e is the radiation efficiency. For SDSS 0836+0054, the most luminous quasar in our sample, it 
takes 20// e-folding times, or ~ 0.8 Gyr, to grow from a 10 Mq stellar black hole, if e = 1. This 
timescale is very close to the age of the universe at z ~ 6 in the A- model. If r/ = 1, it is actually 
longer than the age of the universe in the = 1 model. Note that the black hole could grow much 
faster if the radiation efficiency is very low. Alternatively, the black hole could have formed from 
a much more massive seed black hole, such as a very massive object in Population III (VMOs; 
Bond, Arnett & Carr 1984), or from mergers of smaller black holes (e.g. Kauffmann & Haehnelt 
1999; Menou, Haiman, & Narayanan 2001). 



6.1. Constraint on the Shape of the Quasar Luminosity Function 

The quasar luminosity function at high redshift is a sensitive test of cosmological parameters 
and models of quasar evolution. The observed luminous quasars most probably represent rare 
peaks in the density field at z ~ 6, and hence probe the exponential, high-mass tail of the 
underlying dark matter halo distribution. The slope of the luminosity function is determined by 
both the slope of the halo mass function and the relation between black hole mass (proportional 
to quasar luminosity if the quasar is radiating at the Eddington limit) and dark halo mass. The 
z ~ 6 quasars presented in this paper are extremely rare objects, representing many-u peaks in the 
density field. One would therefore expect their luminosity function to be very steep. We use the 
luminosity distribution of the four quasars in our sample to constrain the slope of the bright end 
of the quasar luminosity function, and compare it with theoretical expectations, assuming that 
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the fluxes of the four observed quasars are not signiflcantly magnified by gravitational lensing. 

Among the four quasars in the complete sample in Table 5, two have M1450 < — 27 (f2 = 1 
model). We use this as a constraint to calculate the expected number of fainter quasars in the 
survey area, assuming a power law luminosity function, ^{L) cx L^, and taking into account the 
survey selection function presented in Figure 8. Wc then compare this expected number with 
the actual number of lower-luminosity quasars observed (two at -M1450 > —27). For /? > —3.9, 
we expect to find less than six quasars at M > —27 in our survey. In this case, the probability 
of observing fewer than three lower luminosity quasars is larger than 5% (assuming Poisson 
statistics). Therefore, /? > —3.9 is consistent with the observed distribution at greater than the 
2-0" level. Note that this constraint encompasses the values of the bright-end slope of the quasar 
luminosity function from Fan et al. (2001b), who found f3 = —2.6 using a sample of quasars at 
2; ~ 4, and Boyle et al. (2000), who found /? = —3.4 at 2; < 2. However, the observed luminosity 
distribution is not consistent with an even steeper luminosity function. The probability drops to 
0.01 for f3 = —4.3, and to 0.001 for /? = —4.6. In other words, a steep quasar luminosity function 
with L oc L~^-^ is ruled out at the ~ 3-a level by the current observations. This constraint simply 
reflects the fact that even after taking into account the lower selection completeness of the survey 
at the faint end, there are not many lower luminosity quasars in our survey at z ~ 6 (only two 
are observed), given that we have discovered two bright ones (M1450 < —27, z* < 19.3). We note 
that this constraint on the slope of the bright end of the luminosity function is quite weak, as it 
is based on only four objects, and depends strongly on our model of the survey selection function 
and the assumption that the observed quasars are not significantly magnified by gravitational 
lensing. In order to obtain a stronger constraint on the slope of the luminosity function, a large, 
and more importantly, a deeper sample of z > 6 quasars is needed. 

Is this constraint on the bright end slope of the quasar luminosity function consistent with a 
simple model of quasar evolution? We can answer this question in the context of specific models 
of quasar formation. Following Haiman & Loeb (1998), we assume that there is a monotonic 
relation between quasar luminosity at high redshift and the mass of the dark matter halo in which 
it resides, and assume further that every black hole shines as a quasar for a given fraction of the 
age of the universe to that redshift (the "lifetime" or "duty cycle" of the quasar). We can then 
match the observed quasar number density at z = 6 to the number density of dark matter halos 
found in the Hubble- volume simulations (Jenkins et al. 2000), if wc assume that the quasars in our 
sample are associated with dark matter halos more massive than Mhaio = 1 - 3 x IO^^Mq. In this 
case, the dark halo mass function at this mass scale has a logarithmic slope between —5 and —6, 
depending on the cosmology and quasar life time (~ 10^ — 10^ years). Note that this minimum 
halo mass is comparable to that estimated using the z ^ relation of Magorrian et al. (1998), 
assuming that the quasars radiate at the Eddington limit. 

These luminous quasars represent about 5-a peaks in the density field for a A-CDM power 

spectrum, and the mass function is extremely steep at this mass. The slope of the dark halo mass 
function is steeper than that of the quasar luminosity function /? > —4.6. This difference in the 



- 21 - 



slopes of the bright end of the quasar luminosity function and the dark halo mass function at the 
corresponding mass scale implies that for high-redshift quasars, the luminosity L cannot simply 
scale linearly with Mhaio! rather, high-mass black holes radiate more efficiently than do lower 
mass ones. Alternatively, the quasar luminosity function can be shallower than the mass function 
if there is large scatter in the L — M relation, or if the fluxes of the observed quasars have been 
significantly amplified by gravitational lensing. Clearly, the above results are oversimplified and 
are based on a number of assumptions regarding both quasar activity and the correspondence 
between quasar luminosity and the mass of the dark matter halo in which it resides. However, 
they illustrate the importance of the measurement of the shape of the quasar luminosity function 
in constraining models of quasar formation and evolution. 



6.2. Quasar Contribution to the Ionizing Background at z 6 

The absence of the Gunn-Peterson trough in SDSS 1306+0356 {z = 5.99) indicates that the 
universe is still highly ionized at z ~ 5.7, presumably by the ionizing photons from quasars and 
star-forming galaxies. The photon emissivity per unit comoving volume required to keep the 
universe ionized at z = 6 in the Q. = 1 model is ( [Madau, Haardt, fc Rees 1999| ), 

*„„(.).10-s-Mpc-3(|)x(i±i)^(2^^^)\ (8) 

where C is the clumping factor of the IGM. Based on the extrapolation of the quasar luminosity 
function from z < 4.5, Madau et al. (1999) conclude that quasars alone do not provide enough 
ionizing photons to keep the universe ionized at z > 5, and that the major contribution to the 
cosmic ionizing background has to be produced by star-forming galaxies (see also pteidel, Pettini 



& Adelberger 2001). Here we calculate the total contribution of the quasar population to the 
ionizing background at z ~ 6, using the quasar density derived in §5.2 and the constraint on the 
bright end slope of the quasar luminosity function in §6.1. 



We calculate the total emissivity per unit comoving volume of ionizing photons from the 
quasar population: 



^Q{z) = „1450 ' (9) 



where e^^^{z) is the total quasar emissivity at 1450A, with units of erg s ^ Hz ^ Mpc ^, assuming 
a luminosity function ^{L, z) 

eQ'°(^)= /^(il450,^)il450f^^l450, (10) 

and n}^^*^ is the total number of ionizing photons for a source with a luminosity of 1 erg s^^ Hz~^ 
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at 1450A. We assume the quasar SED following Madau, Haardt, &; Rees (1999). 

L(i/) oc < 



2500A< A < 4400A; 

1050A< A < 2500A; (11) 
z/-^-^ A < 1050A, 



and integrate over the energy range from 1 to 4 Rydberg, with a cutoff at 4 Rydberg because of 
Hell absorption. 

We assume that the luminosity function of quasars at z ~ 6 is described by a double power 
law, and write it as a function of absolute magnitude M1450: 

*(Mi45o) = ^q0.4(/3i+1)(Mi4B0-M*45J ^ iQ0A{/32+l){Mii50-M*^^^) ' (^2) 

where M^^^^q is the characteristic luminosity. 

None of the current quasar surveys at z > 4 constrains the faint end slope of the quasar 
luminosity function. For example, Fan et al. (2001b) only compute the luminosity function at 
-^"1450 < —25.5. Therefore, the characteristic luminosity M^^r^Q and the faint end slope /?2 are both 
unknown at high redshift, and we simply assume the faint end slope of the luminosity function to 
be the same as determined by Boyle et al. (2000) at z < 2, P2 = —1.58. Thus, for a given bright 
end slope Pi, we can derive ^* by requiring JZ^'^ ^{M)dM = 1.1 x 10~^ Mpc~^ , as determined 
in §5.2. 

In §6.1, we demonstrated that the bright end slope Pi of the quasar luminosity function 
at z ~ 6 is consistent with that measured at lower redshift (Pi = —2.58, Fan ct al. 2001b, and 
Pi = —3.43, Boyle et al. 2000), and that a much steeper slope can be ruled out {Pi < —4.3 is ruled 
out at the 99% confidence level, and Pi < —4.6 with 99.9% confidence). In Figure 10, we calculate 
the quasar ionizing photon density integrated over the full luminosity function as a function of 
M^450 for a series of values of Pi, assuming P2 = —1.58, and normalizing to the quasar density at 
-^1450 < —26.8. The quasar ionizing photon density is compared with what is required to keep the 
universe ionized at z ~ 6 (heavy solid line). 

From Figure 10, it is evident that with Pi > —3.5, quasars would not produce enough photons 
to keep the universe ionized at z 6, unless M^^^q > —23, i.e., only if the ionizing background is 
dominated by low luminosity AGNs can the quasar luminosity function provide enough ionizing 
photons. This would imply a large negative luminosity evolution of the luminosity function. The 
universe could be ionized with luminous quasars only if the bright end slope of the luminosity 
function were much steeper than that measured at low redshift; but such a steeper slope is 
not consistent with our sample. This calculation underlines the importance of searching for 
faint quasars and studying the shape of the high-redshift quasar luminosity function at the low 
luminosity end. While it is unlikely that quasars provide the majority of photons that ionized the 
universe at z ^6, it still remains an open question. 
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7. Summary 

In this paper, we present three new quasars at 2; > 5.8, selected from ~ 1550 deg^ of multicolor 
imaging data from the Sloan Digital Sky Survey. The quasar candidates axe selected as i-dropout 
objects, with red i* — z* colors. They are separated from red cool dwarfs using follow-up J-band 
photometry. The three new quasars, with redshifts of 5.82, 5.99 and 6.28, are the highest-redshift 
objects with spectral confirmation known to date. 

The spectra of the new quasars show very strong Lya absorption blueward of the Lya 
emission, with more than 90% of the flux absorbed. The amount of absorption increases rapidly 
toward Zabs ~ 6. For the z = 6.28 quasar SDSS 1030+0524, the flux is consistent with zero in 
a region of 300A immediately blueward of Lya emission. It suggests a tentative detection of a 
complete Gunn-Peterson trough, indicating that at 2; ~ 6 the universe is close to the reionization 
epoch. 

The discovery spectra show tentative detection of the NV emission line, suggesting the 
existence of large amount of heavy metals in the gas around the quasars. A J-band spectrum of 
SDSS 1030+0524 shows strong CIV emission. The limits on the NV/CIV and NV/Hell line ratios 
are consistent with the gas around the quasar having super-solar metallicity. 

The three new quasars described in this paper, plus the z = 5.8 quasar SDSS 1044-0125 (Fan 
et al. 2000c), form a complete color-selected, flux- limited sample at z ~ 6. We calculate the color 
selection completeness of this sample, and derive the spatial density of luminous quasars at 2; 6, 
p(Afi45o < -26.8) = 1.1 X 10-9 Mpc-3, for O = 1 and Hq = 50 km s'^ Mpc'^. This density is 
about a factor of two lower than that at 2; ~ 5, and is consistent with an extrapolation of the 
observed redshift evolution of quasars at 3 < 2; < 5. 

Using the luminosity distribution of the sample, we show that the luminosity function is 
shallower than *(-L) oc L~^-^ at the bright end. A larger and deeper sample is needed to better 
constrain the shape of the quasar luminosity function. We confirm that the quasar population is 
unlikely to provide enough photons to ionize the universe at 2; ~ 6, unless the luminosity function 
is much steeper than P2 = —1.58 at the faint end. 

The black hole masses of these quasars are probably several times lO'^M©. The quasars are 
likely to reside in very massive systems, with the minimum mass of host dark halos ~ IO^'^Mq. 
These massive dark halos represent rare peaks in the density field at high redshift, and arc in the 
steep tail of the mass function, with a slope that is appreciably steeper than that of the quasar 
luminosity function. 

Follow-up observations of the luminous high-redshift quasars described in this paper will 
provide excellent probes of galaxy formation and IGM evolution at high redshift. In Paper II, we 
present high-resolution, high-S/N optical spectroscopy of these quasars. These observations are 
compared with detailed cosmological simulations to constrain models of reionization in Fan et al. 
(2001c). X-ray observations (e.g., Brandt et al. 2001) will reveal the status of the environment 
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very close to the quasar central engine; follow-up IR spectroscopy can provide more diagnostics 
on the metallicity of the quasars and reveal the possible existence of strong associated absorption 
(e.g., [Maiolino et al. 2001 ); and detections of the sources at sub-millimeter wavelengths will shed 



light on the possible connections between quasar activity and starbursts (e.g., Carilli et al. 2001). 

The total area of the SDSS survey is 10,000 deg^. Assuming the quasar luminosity function 
of Fan et al. (2001b) and normalizing it to have the spatial density at z ~ 6 found in this paper, 
we expect to detect one z ~ 6.6 quasar with z* ^ 20 {Mi^^q ~ —28) in the entire survey. A quasar 
at this redshift has i* — z* and z* — J colors similar to those of early T dwarfs. This is the highest 
redshift of objects we expect to find using solely SDSS imaging data. 

The Sloan Digital Sky Survey (SDSS) is a joint project of The University of Chicago, 
Fermilab, the Institute for Advanced Study, the Japan Participation Group, The Johns Hopkins 
University, the Max-Planck-Institute for Astronomy (MPIA), the Max-Planck- Institute for 
Astrophysics (MPA), New Mexico State University, Princeton University, the United States Naval 
Observatory, and the University of Washington. Apache Point Observatory, site of the SDSS 
telescopes, is operated by the Astrophysical Research Consortium (ARC). Funding for the project 
has been provided by the Alfred P. Sloan Foundation, the SDSS member institutions, the National 
Aeronautics and Space Administration, the National Science Foundation, the U.S. Department 
of Energy, the Japanese Monbukagakusho, and the Max Planck Society. The SDSS Web site is 



http:/ /www. sdss.org . This publication makes use of data products from the Two Micron All Sky 
Survey, which is a joint project of the University of Massachusetts and the Infrared Processing 
and Analysis Center/California Institute of Technology, funded by NASA and NSF. We thank the 
staff at APO, Keck, Calar Alto, and UKIRT, for their expert help, Gemini staff scientist Marianne 
Takamiya for carrying out the Germini/NIRSPEC service observation. We thank Ed Turner, 
Peng Oh, Julian Krolik, Kristian Finlator, John Bahcall and an anonymous referee for helpful 
discussions and comments. We acknowledge support from NSF grant PHYOO- 70928 and a Frank 
and Peggy TapUn Fellowship (XF), NSF grant AST-0071091 (MAS), and NSF grant AST-9900703 
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Table 1. Summary of Photometric Runs 



run 


strip 






( FWHM(i) ) (") 


( FWHM(2;) ) ( ) 


745 


10 N 


51.1 


65.0 


1.23 




1.25 


752 


10 S 


-17.6 


51.0 


1.43 




1.50 


756 


10 N 


—63.1 


51.0 


1.34 




1.38 


1140 


9 N 


—30.5 


17.1 


1.44 




1.47 


1231 


9 S 


-14.2 


14.7 


1.51 




1.53 


1239 


10 S 


-63.1 


-39.7 


1.41 




1.46 


1241 


9 S 


—45.0 


—38.9 


1.49 




1.57 


1241 


9 S 


-24.8 


-19.0 


1.74 




1.76 


1331 


36 S 


—56.4 


—36.0 


1.18 




1.18 


1331 


36 S 


—30.2 


—18.0 


1.38 




1.36 


1332 


36 S 


-18.1 


-10.2 


1.52 




1.46 


1336 


42 N 


23.4 


36.1 


1.51 




1.50 


1339 


42 S 


23.3 


36.0 


1.41 




1.38 


1345 


36 N 


-52.6 


-8.0 


1.42 




1.42 


1345 


36 N 


1.3 


73.9 


1.23 




1.26 


1350 


37 S 


-53.7 


14.5 


1.54 




1.54 


1356 


43 N 


22.8 


36.5 


1.64 




1.62 


1359 


43 S 


27.7 


37.2 


1.61 




1.57 


1402 


37 N 


-51.0 


-34.1 


1.24 




1.22 


1412 


37 N 


—23.6 


33.2 


1.28 




1.28 


1449 


1 N 


2.9 


13.4 


1.46 




1.51 


1450 


37 N 


-34.2 


—23.3 


1.53 




1.54 


1462 


11 S 


-56.1 


40.6 


1.27 




1.35 


1468 


32 N 


10.4 


29.1 


1.37 




1.35 


1469 


32 S 


7.5 


29.4 


1.49 




1.48 


1478 


12 S 


25.6 


55.8 


1.43 




1.49 


2074 


35 S 


—34.6 


—25.2 


1.51 




1.56 


2076 


35 N 


—56.2 


—33.9 


1.33 




1.35 


2078 


36 S 


-9.4 


8.5 


1.47 




1.49 


2125 


12 S 


-64.5 


-17.1 


1.53 




1.58 


2126 


12 N 


-64.2 


0.2 


1.30 




1.36 


2131 


34 S 


-60.3 


-42.6 


1.41 




1.50 


2134 


36 S 


8.4 


23.5 


1.48 




1.41 


2137 


34 N 


-59.3 


-19.6 


1.37 




1.43 


2138 


34 N 


-22.7 


-9.5 


1.51 




1.55 


2140 


34 S 


-45.4 


-35.3 


1.69 




1.76 


2141 


10 S 


-38.0 


-17.2 


1.42 




1.51 


2190 


12 N 


2.3 


28.2 


1.47 




1.54 


2190 


12 N 


32.9 


51.7 


1.48 




1.53 
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Tablc 2. Summary of Follow-up Results 



number of objects percentage 



z > 5.8 quasars 


4 


3.3% 


T dwarfs 


11 


9.1% 


M/L dwarfs by spectroscopy 


10 


8.3% 


M/L dwarfs by photometry 


60 


49.6 % 


BAL quasar 


1 


0.8% 


false detection 


35 


28.9% 



TOTAL 121 100% 



Table 3. Photometric Properties of Three New z > 5.8 Quasars 



object 


redshift 


i* 


z* 


J 


SDSS run 


J083643.85+005453.3 
J103027.10+052455.0 
J130608.26+035626.3 


5.82 ± 0.02 
6.28 ± 0.03 
5.99 ± 0.03 


21.04 ± 0.08 

23.23 ± 0.43 
22.58 ± 0.26 


18.74 ± 0.05 
20.05 ± 0.10 
19.47 ± 0.05 


17.89 ± 0.05 
18.87 ± 0.10 
18.77 ± 0.10 


1239 
2125 
2190 



The SDSS photometry {i*,z*) is reported in terms of asinh magnitudes on the AB system. The 
asinh magnitude system is defined by Lupton, Gunn & Szalay (1999); it becomes a linear scale in 
flux when the absolute value of the signal-to-noise ratio is less than about 5. In this system, zero 
flux corresponds to 24.4 and 22.8, in i*, and z*, respectively; larger magnitudes refer to negative 
flux values. The J magnitude is on the Vega-based system. 
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Table 4. Average Lya Absorption oi z > 5.8 Quasars 



object 


redshift 


Da 


Db 


^abs 


Transmitted flux ratio 














J104433.04-012502.2 


5.80 


0.91 ± 0.02 


0.95 ± 0.02 


5.5 


0.088 ± 0.004 


J083643.85+005453.3 


5.82 


0.90 ± 0.02 


0.91 ± 0.02 


5.5 


0.099 ± 0.014 


J130608.26+035626.3 


5.99 


0.92 ± 0.02 


0.95 ± 0.02 


5.7 


0.069 ± 0.014 


J103027. 10+052455.0 


6.28 


0.93 ± 0.02 


0.99 ± 0.01 


6.0 


0.003 ± 0.020 



Table 5. Continuum Properties oi z > 5.8 Quasars in the Complete Sample 



object 


redshift 


?TZ1280 




Mi45o (O-model) 


Mi45o (A-model) 


E{B-l 


J083643.85+005453.3 


5.82 


18.88 


18.81 


-27.62 


-27.88 


0.050 


J103027.10+052455.0 


6.28 


19.73 


19.66 


-26.89 


-27.15 


0.023 


J104433.04-012502.2 


5.80 


19.28 


19.21 


-27.15 


-27.50 


0.054 


J130608.26+035626.3 


5.99 


19.61 


19.55 


-26.93 


-27.19 


0.028 
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Figure 1. i* — z* vs. z* color-magnitude diagram for the i-dropout sample. Colors and 
magnitudes are plotted as asinh magnitudes measured by SDSS imaging. The symbols represent 
the classification in Table 2: circles: T dwarfs; triangles: L dwarfs; stars: z > 5.8 quasars; square: 
BAL quasar. Filled symbols are objects with S/N in the i band higher than 4; open symbols are 
objects not detected in the i band at the 4-a level. The median track of simulated i* — z* color 
and z* magnitude for quasars with M1450 = —27 is also shown as a function of redshift, with plus 
signs every 0.1 in redshift. For comparison, the data for a random sample of 50,000 high-latitude 
stars are also shown as dots. The dashed line shows the cut i* — z* > 2.2 that we use to select 
high-redshift quasar candidates. 
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Figure 2. i* — z* vs. z* — J color-color diagram for the i-dropout sample. The symbols are the same 
as in Figure 1. The median track of simulated quasar colors is shown as a function of redshift. The 
survey selection criteria are illustrated by the shaded area. For comparison, colors of SDSS-2MASS 
stars in a 50 deg^ area at high latitude are also shown. The arrow at the lower-right corner indicates 
the reddening vector for quasars at z = 6 with E{B — V)= 0.1. 




Figure 3. Projection of the survey area in J2000 equatorial coordinates. Note that for some SDSS 
stripes, only one of the two scans is covered by the current survey. The stripes are each 2.5° wide. 
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Figure 4. z finding charts of tlie tliree new z > 5.8 quasars. For SDSS 0836+0054, the SDSS z 
image is shown. For SDSS 1030+0524 and SDSS 1306+0356, the ARC 3.5m z images taken with 
SPICAM are shown (180 seconds exposure). The size of the finding chart is 160". North is up and 
East is left. 
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Figure 5. The discovery spectra of the three new z > 5.8 quasars, taken with the ARC 3.5m 
telescope and DIS spectrograph. The exposure time is 3600 seconds for each quasar and the 
spectral resolution is about 20A. The spectra have been normalized to the photometry in the z 
filter. 
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Figure 6. Combined optical + near IR spectrum of SDSS 1030+0524. The optical spectrum is a 
3600 second exposure taken with ARC 3.5m telescope (same as in Figure 5). The near IR spectrum 
is a 3000 second exposure taken with Keck/NIRSPEC. The resolution of the near IR spectrum is 
R ~ 1500. 
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Figure 7. Evolution of the transmitted flux ratio in the Lya forest region, averaged over the 
wavelength range (1 + Zabs — 0.1) x 1216A ^ A ^ (1 + Zahs + 0.1) x 1216A. The vertical error bars 
represent the photon noise. Note that for SDSS 1030+0524, the flux is consistent with zero in the 
Lya region immediately blueward of the Lya emission. 




Figure 8. The selection probability oi z > 5.8 quasars as a function of redsliift and luminosity, for 
the = 1 and the A models. The probability contours of 5%, 20%, 40%, 60% and 80% are shown. 
The large dots represent the locations of the 4 quasars in the sample. 
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Figure 9. The evolution of quasar comoving spatial density at M1450 < —26.8 in the = 1 model. 
The large dot represents the result from this survey. The dashed and dotted lines are the best-fit 
models from Fan et al. (2001b) and Schmidt et al. (1995), respectively. The solid line is the best-fit 
model from the 2dF survey (Boyle et al. 2000) at 2; < 2.5. 
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Figure 10. Comoving emission density rate of hydrogen Lyman continuum photons from quasars 
integrated over the full luminosity function, compared with the minimum rate (heavy solid line) 
that is needed to fully ionize a universe at z = 6, in the $1 = 1 model, with a clumping factor 
of 30, and the baryon fraction Q^h^ = 0.02. The recombination time is assumed to be faster 
than the Hubble time. The quasar emission rate is calculated for different bright end slopes of 
the quasar luminosity function and turnover luminosity Mf^^Q, assuming the same faint end 
slope P2 = —1-58 and the total density of luminous quasars matched to that found in this paper, 
p{z = 6,Mi45o < -26.8) = 1.14 x lO'^ Mpc-^. 



